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The adsorption of benzene and C60 on graphene and boron nitride (BN) is studied using density 
functional theory with the non-local correlation functional vdW-DF. By comparing these systems 
we can systematically investigate their adsorption nature and differences between the two functional 
versions vdW-DFl and vdW-DF2. The bigger size of the C60 molecule makes it bind stronger to 
the surface than benzene, yet the interface between the molecules and the sheets are similar in 
nature. The binding separation is more sensitive to the exchange variant used in vdW-DF than to 
the correlation version. This result is related to the exchange and correlation components of the 
potential energy curve (PEC). We show that a moderate dipole forms for C60 on graphene, unlike 
for the other adsorption systems. We find that the corrugation is very sensitive to the variant or 
version of vdW-DF used, in particular the exchange. Further, we show that this sensitivity arise 
indirectly through the shift in binding separation caused by changing vdW-DF variant. Based on our 
results, we suggest a concerted theory-experiment approach to assess the exchange and correlation 
contributions to physisorption. Using DFT calculations, the corrugation can be linked to the optimal 
separation, allowing us to extract the exchange-correlation part of the adsorption energy. Molecules 
with same interfaces to the surface, but different geometries, can in turn cast light on the role of 
van der Waals forces. 

PACS numbers: 68.43.-h,71.15.Mb,73.90.+f,81.05.ue 



I. INTRODUCTION 

The development of methods that include van der 
Waals forces in density functional theory (DFT) has 
made it possible to describe an abundance of impor- 
tant material systems — like water j 1 * 2 - DNA^ molecular 
crystals^— and metal-organic frameworks^ — at the 
electronic level even without use of empirical parame- 
ters. Such materials can be categorized as sparse matter, 
because of the essential role of the low (electron) density 
regions or "voids" in these materials. These regions make 
it paramount to go beyond semi-local approximations like 
the generalized-gradient approximation (GGA)^£r— and 
properly account for van der Waals interactions. A good 
description of sparse matter is important for developing 
new nanotechnology. The softer van der Waals bonds 
makes it easy for fragments to move, like in the self as- 
sembly of molecules on surface a 17 i 18 or for moving pieces 
of nano-mechanical devices^ 

Sparse matter constitutes a research-frontier for DFT 
and DFT-based modelling. 20,21 Practical and robust 
schemes to describe sparse matter are all fairly recently 
developed and several competing and complimentary 
methods exist. Two prominent approaches are vdW- 
corrected DFT and the use of non-local correlation func- 
tionals. The first approach adds a semi-empirical pair- 
potential correction accounting for van der Waals forces 
on top of semi-local DFTi22r— There exist many vari- 
ation thereof and refinements, like adjusting the pair- 
potentials by using charge-density input. 26 Non-local cor- 
relation functionals aim to stay fully within DFT keep- 
ing the density as the key variable. The van der Waals 
density functional (vdW-DF) metho d 21 ! 27 " — for approxi- 



mating the exchange-correlation energy has led to a series 
of such functionals that lack empirical input. The vdW- 
DFl^ introduced in 2004, and vdW-DF2,22 introduced 
in 2012, emphasise a constraint-based design of the den- 
sity fluctuation propagator while the precursor for 
layered system a 2 : 28 emphasise an anisotropic screening 
account. The vdW-DFs vanish seamlessly in the uniform 
limn— and exhibits the correct 1/r 6 asymptotic limit 
between molecular dimers^ The prediction of various 
sparse matter methods can differ quite much. 

Hand in hand with the development and refinements of 
sparse matter methods, it is essential to identify exper- 
iments and quantum-chemistry calculations that serves 
to test and elucidate these methods for various physical 
effects and across length scales£i2£r— One approach is 
to construct a representative data set that account for 
many physical effects; 41 - - — like the S22 data set42 or the 
more recent and larger S66. 44 The average performance 
of methods can then be compared to each other. Another 
approach is to identify particular systems and associated 
experiments that reveal several properties of sparse mat- 
ter binding^ An example is to calculate potential en- 
ergy curves (PEC) of H 2 on Cu(lll) and compare with 
the experimental ones constructed from backscattering 
experiments! 46 ' 47 Since the balance between the differ- 
ent contributions to PEC change with separation, much 
information can be extracted. 

In this paper, we show that the details of adsorption, 
especially the predicted corrugation for adsorbate dy- 
namics is a discriminator for assessing the balance be- 
tween exchange and correlation components in a non- 
local functional design. We substantiate our conclusion 
by comparing the results of calculations using vdW-DFl 
and vdW-DF2 and associated exchange choices for a 
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FIG. 1: Adsorption systems studied. Upper left: benzene on 
graphene on top site. Upper right: benzene on BN on top- 
N site. Lower left: Repeated supercells of C60 on graphene. 
Lower right: C60 on BN on top-N-rot site (viewed from below 
BN sheet). 



group of adsorbate systems that are closely related, yet 
the nature of the binding differs in a well-defined manner. 

We study the adsorption of benzene and C60 on 
graphene and BN. We show how these systems allow us 
to span, in a systematic manner, a wide range of different 
contributions to sparse matter binding. First, the inter- 
face between a benzene and graphene (or BN) is very sim- 
ilar to that of C60 on graphene (or BN), both in terms 
of contact area and character of the local bonds. Yet 
C60 is bigger and the difference in adsorption is mostly 
induced by the increased van der Waals attraction of 
the C60 molecule. Second, boron nitride and graphene 
are both flat one-atom layer thick sheets with eight va- 
lence electrons per unit cell, but whereas graphene is a 
semi-metal, BN is an insulator. The BN sheet also has 
a somewhat different charge-density landscape than the 
graphene sheet. 

We show that the corrugation is very sensitive to the 
binding separation and therefore very sensitive to the ver- 
sion of correlation and exchange in particular used in 
vdW-DF. Yet at a given separation to the surface, the 
corrugation is largely independent of the functional ver- 
sion used. The shorter binding separation for C60 adsor- 
bants compared to benzene increase corrugation in the 
same indirect manner. 

Based on our results and analysis, we suggest a com- 
bined theory-experiment approach that can lead to in- 
sight into the nature of sparse matter binding and 




FIG. 2: The vdW-DF kernel. 
Dion et.al.— 



Figure similar to that of 



the magnitude of the different terms in the exchange- 
correlation functional. 

The plan of the paper is as follows: After the intro- 
duction we summarize vdW-DF method and functional 
versions and variants and outline computational details 
specific to this study. Section [Hi] presents our results for 
adsorption of benzene and C60 on graphene and boron 
nitride. In section IIVI we analyze the various compo- 
nents to this binding, both the importance of GGA ex- 
change choice and vdW-DF correlation version. We also 
analyze the role of different density separations and den- 
sity regimes in the non-local correlation of vdW-DF. Sec- 
tion [V] compares the corrugation of the different systems 
and for the different functionals. Section IVII discuss the 
possibility of using the selected systems to assess the 
quality of non-local correlation functionals and sparse- 
matter methods in general. Finally, we summarize our 
conclusions. 



II. VDW-DF CALCULATIONS 

The vdW-DF method is designed to handle both short- 
range covalent bonds and long-range van der Waals in- 
teractions without introducing damping functions or em- 
piricism. Its exchange-correlation functional combines 
exchange at the GGA level, ££ dwDF = E^ GA with the 
combination of correlation in the local density approx- 
imation (LDA) and non-local correlation: £* dWDF = 
non-local correlation is responsible for 
the van der Waals interaction. 

The non-local correlation term vanishes in the uniform 
limit and takes the form of a six dimensional integral 

Ef[ n ] = ^Jd 3 rJ d 3 r'n(r)0[n](r,r>(r'). (1) 

The kernel can be expressed in terms of a universal kernel 
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</>(r, r') = 4>\{d,d') that depends on two dimensionless 
length scales d — go( r i) r i2 an d d' = qa(r 2 )ri 2 where 
1*12 = r i — r2 1 and go is a modulation of the local Fermi 
vector kp(n) that account for the local response. By 
depending on gradients as swell as local density, this local 
response reflects a broader dependence on the density 
variation through a plasmon-pole approximation/ 9 - 

The universal kernel <pi kernel can also be tabulated in 
terms of D = (d+d')/2 and S = (d- d')/(d + d'). Fig. [2] 
shows the kernel as a function of D for three selected 8. 
The 6 parameter is a measure of how different the local 
response of the two density regions are. 

The GGA exchange energy can be expressed as a mod- 
ulation F x of the LDA exchange energy e^ DA (rt(r)), as 
follows 

£ X GGA = J d 3 rn(r) e ^A (?l(r))Fx(s) (2) 

Here s = (|Vn|/2n)/fcF is a measure of how rapidly the 
density is changing compared to length scale set by the 
local Fermi wave vector. In the small s limit, the modu- 
lations of the GGAs go like F x (s) — 1 ~ s 2 , while in the 
large s limit they differ significantly both quantitatively 
and qualitatively. 

A. Functional versions 

Because the nature of the adsorption systems consid- 
ered here differ in a systematic manner, it becomes inter- 
esting to study how the predictions of vdW-DF changes 
with correlation version and choice of exchange partner. 
For vdW-DFl correlation^ 9 , (the version of 2004), we will 
in, in addition to the canonical choice of revPBE, use 
C094 8 . PW86r^i2 and optPBE^ . as exchange partner. 
For vdW-DF2 correlation,— we use C09, in addition to 
the canonical choice of PW86r. 

Standard vdW-DFl with the revPBE 51 exchange func- 
tional as companion to its non-local correlation have 
shown an impressive robustness, giving good results 
for binding between molecules of various sizes and 
properties^— ^ i 21 ' 29 i 31 ' 52 " — adsorption onto insulating 
surfaces/^— graphene/ 9 - - — and metals ^ i 46 i 47 i 62- - 65 anc j 
between various layered compounds . 66 ' 67 Yet, its accu- 
racy for sparse matter does not provide the full accu- 
racy of traditional DFT calculations for typical covalcntly 
bonded syste, where the generalized-gradient approxima- 
tion excels. Some of this inaccuracy can be attributed to 
the overly repulsive revPBE™ exchange functional caus- 
ing overestimated separations. This exchange functional 
was chosen because many other standard choices like 
PBE 15 tend to produce spurious binding arising from ex- 
change effects^^iS 9 . 

Some sparse matter systems are also particularly chal- 
lenging, like those with weak chemisorption. These sys- 
tems are highly sensitive to the balance between the dif- 
ferent terms in the exchange-correlation functional/ -/ 6 - 
The overestimation of binding separations in vdW-DF 



can make it miss a van der Waals-induced charge trans- 
fer, a signature of weak chemisorption. 

Refinements of the vdW-DFl functional have been sug- 
gested, in particular it has been suggested to replace its 
exchange functional with a less repulsive one. Cooper 48 
designed an exchange functional C09 based on differ- 
ent fundamental constraints than revPBE, following the 
principles that underpins the introduction of the PBEsol 
functional/- C09 is less repulsive than revPBE, but be- 
comes similar in the large s-limit, thereby minimizing 
contributions to unphyscial exchange binding. 

Klimes and coworkers^ - reparameterized a set of ex- 
change functionals by optimizing their performance to 
the S22 data set 42 One of these, the optPBE functional, 
corresponds to a reparameterization of PBE and is there- 
fore somewhat similar to revPBE. This functional will be 
included in our study. 

Recently, a new version of vdW-DF called vdW-DF2 
has been developed/ 2 - This functional sets the <?o( r ) pa- 
rameter by using a different description of the plasmon 
that is suggested by an analysis of the exchange- variation 
in the high-density limit £LM. The updated plasmon de- 
scription makes vdW-DF2 less sensitive to low-density 
regions and generally reduces the non-local correlation 
energy. In addition to adjusting the version of corre- 
lation, this functional updates the account of exchange 
by relying on a refitted version of PW86, here labeled 
PW86ri 12 i 49 This exchange functional agree well with the 
exchange of Hartre-Fock calculations and avoids spurious 
exchange bindingi 49 i 68 i 81 This functional is also less re- 
pulsive than revPBE at binding separations. 

In almost every case, the modifications of vdW-DF de- 
scription reduces binding separations compared to the 
results of vdW-DFl. For the S22 data set, 42 vdW-DFl 
with C09 exchange and vdW-DF2 produces good results; 
unsurprisingly so does using vdW-DFl with optPBE, 
since it is was fitted to this data set. When compar- 
ing the versions for a broader class of systems, we find 
that no version clearly outperforms the others. 

To exemplify this observation, vdW-DF2 stands out 
by predicting a potential energy curve for H2 on Cu(lll) 
that compares well with that generated from experimen- 
tal data ; 46 ' 47 both in terms of optimal separation and ad- 
sorption energy. However it underestimates the asymp- 
totic form. For C60 crystals/ vdW-DFl predicts a cohe- 
sive energy 1.6 eV in good agreement with experimental 
measured interval 1.6-1.9. vdW-DFl using C09 exchange 
functional predicts an excellent lattice constant and a co- 
hesive energy of about 2 eV. In contrast vdW-DF2 pre- 
dicts a cohesive energy of merely 1.3 eV. For adsorp- 
tion of C60 on Au(lll), the combination of vdW-DF2 
and C09 has been found to yield the best results/ 1 - For 
this system, neither vdW-DFl nor vdW-DF2 are able to 
predict the right balance between repulsive and attrac- 
tive terms that results in the experimental finding of a 
weak chemisorption. 70 We find that vdW-DFl predicts 
the best binding energies for both benzene and C60 on 
graphene. 
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Bjorkman and coworkers recently asked the question 
"are we van der Waals ready?" in a large survey of the 
ability of electronic structure methods, including vdW- 
DF1 and vdW-DF2, to accurately describe binding of 
layered systems 4i None of the DFT functional meth- 
ods they tested are able to consistently predict accurate 
binding energies, separations, and elastic constants. For 
instance both vdW-DFl and vdW-DF2 in general over- 
estimates interlayer separations. 

To get a better understanding of why different func- 
tionals prevail for different systems, a more systematical 
analysis of their binding mechanisms are needed. A sug- 
gestion to test performance by comparing with exper- 
iments that reveal many interaction properties for the 
same one system was launched in Refs. l4q|47l Here, we 
seek to map out how both the exchange and correlation 
functionals contribute to PECs for related molecules of 
different sizes and based on this analysis propose an addi- 
tional adsorption-based strategy to refine functional test- 
ing. Some of the analysis approaches taken here can be 
extended to other van der Waals bonded systems, or di- 
rectly applied to other non-local correlation functionals 
like those of Vydrov and Voorhis^r— 



B. Computational method 

The vdW-DF adsorption study presented here relies on 
a post-GGA procedure similar to many previous vdW- 
DF studies: the charge density is first generated in a 
self-consistent PBE calculation with ultrasoft pseudopo- 
tentials using the Dacapo software package^ next, in 
a post-processing procedure, we calculate the semi-local 
exchange at the GGA level and non-local vdW-DF cor- 
relation. These terms replace those of the self-consistent 
calculation. The non-local correlation part is calculated 
with an in-house real-space code presented and discussed 
in Ref . [D- 

This study is similar to previous vdW-DF adsorption 
studies detailed in Refs. I 8.60.83 and resolves computa- 
tional issues like grid sensitivity and handling of super- 
cells in the same manner. Here, the plane- wave cutoff 
is set to 500 eV and a k-sampling of 2 x 2 x 1 is used. 
The supercell consists of 6 x 6 graphene (boron nitride) 
primitive unit cells. This size makes inter-adsorbate in- 
teraction negligible even for the case of C60 adsorbates 
(see lower left of Fig. [IJ , since the non-local interaction 
between adsorbates is explicitly corrected for. 

It is useful to consider the part of the total energy in 
vdW-DF that contains one-particle kinetic, electrostatic, 
local correlation, and exchange Eq separately from the 
remaining non-local correlation E~. In a purely van der 
Waals bonded system, the former part typically gives 
rise to a repulsive contribution to the potential energy 
curve (PEC), while the non-local correlation part, an at- 
tractive. In calculating the contributions to the PEC 
between a molecule and surface, we here calculate the 
former by subtracting off the energy off a reference sys- 



tem where the molecule is far from the surface 

AE = E (h)-E (h oo ); (3) 

in practice hoc = 9 A is sufficient. While for the 
contribution arising from non-local correlation, which 
have longer range, we subtract off the energy of isolated 
fragmentsiSl 

A^c' = ^c.fullC 1 ) — ^c 'sheet — ^"c 'molecule ■ (4) 

III. ADSORPTION RESULTS 

We describe the adsorption results of our vdW-DF 
study of adsorption of benzene and C60 on graphene and 
boron nitride. We present the binding energies and sheet- 
to-molecule separation for the different systems and for 
the different functionals versions. We also calculate the 
dipole induced by adsorption and show that only C60- 
on-graphene induces a significant dipole at binding sep- 
arations. 



A. Binding energy and separation 

To determine the optimal site and corresponding bind- 
ing energy and separation, we optimize the separation 
to the sheet for each high-symmetry site. Fig. [3] shows 
the six high symmetry sites considered for benzene on 
graphene. The molecule is assumed to lay flat on the 
sheet. On boron nitride, we consider eight high symme- 
try sites, since the center ring can be on-top of a nitrogen 
(top-N/top-N-rot) and on-top of a boron {top-B/top-B- 
rot). For C60 we consider both configurations where the 
carbon hexagon or pentagon face the sheet and find that 
the hexagon is the preferred one. For benzene, the opti- 
mal site is the top on graphene and top-N on boron ni- 
tride. The optimal hexagonal configuration are rotated 
30° degrees about the surface normal z for C60 compared 
to the optimal one for benzene (top-rot/top- N-rot) . The 
same optimal sites were predicted by the considered ver- 
sions and variants of vdW-DF. 

Figure |4] shows PECs for benzene and C60 on graphene 
as calculated in vdW-DFl and vdW-DF2. Table U dis- 
plays the adsorption energies and the optimal molecule- 
to-sheet separations for the optimal sites including ad- 
sorption on BN and use of non-canonical exchange part- 
ners. 

The benzene adsorption energies on graphene almost 
match those on BN, while C60 binds slightly stronger 
to graphene. The optimal benzene-on-graphene sepa- 
ration is slightly larger than the benzene-on-BN, while 
the C60-on-graphene and C60-on-BN separations are the 
same. That benzene and C60 does not follow identical 
trends could stem from a moderate charge transfer be- 
tween graphene and C60, compared to a much smaller 
one for the C60 on BN. The difference in charge transfer 
is discussed in the next subsection. 




FIG. 3: Adsorption sites considered for benzene on graphene. 
The -rot sites are rotated 30°. The top site is the optimal 
one. 
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FIG. 4: The potential energy curves (PEC) for benzene (yel- 
low curves) and C60 (black) for adsorption on graphene as 
obtained with vdW-DFl (full) and vdW-DF2 (dashed). 



vdW-DFl predicts adsorption energies very close to 
the experimental for both benzene and C60 on graphene. 
vdW-DF2 somewhat underestimates the experimental 
value. vdW-DFl correlation combined with the non- 
canonical exchange variants considered here all result in 
overestimated adsorption energies, while vdW-DF2 with 
C09 exchange underestimates the adsorption energies. 

Both vdW-DFl and vdW-DF2 describe well the ex- 
pected increased van der Waals attraction for the larger 
system. vdW-DFl predicts that C60 binds 73% stronger 
than benzene on graphene, vdW-DF2 predicts 67%, and 
the experimental number is 70%. Note though, the spe- 
cific factor change if a different exchange functional is 
used. 

We are unaware of explicit experimental adsorption en- 
ergies for benzene and C60 on BN. Reinke and coworkers 
did find that those C60 molecules not associated with 
defects desorbed at slightly lower temperature on boron 
nitride than on graphene^ This observation is consis- 
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TABLE I: Adsorption energy i? a ds and separation h between 
the molecule and graphene (closest atom) as obtained for dif- 
ferent variants of vdW-DF. Energies are compared to the ex- 
perimental for adsorption ongraphene. 



benzene on graphene benzene on BN 



Functional 


h [A] i?/ ads [eVJ 


h [A] £ ads 


eV] 


VUVV 171 A 


q r n aq 

o.u u.^ty 


3.5 0.48 




vdW-DFl(C09) 


3.25 0.59 


3.2 0.59 




vdW-DFl(optPBE) 


3.4 59 


3.35 0.58 




vdW-DFl(PW86r) 


3.35 0.66 


3.35 0.66 




vdW-DF2 

V V L ti J J J_ £d 


3.5 0.43 


3.4 0.43 




vdW-DF2fC09') 


3.4 0.34 


3.3 0.34 




JJjXp. 


U.OU3ZU.UO 








C60 on graphene 


C60 on BN 


Functional 


h [A] £ ads [eV] 


h [A] £ ads 


eV] 


vdW-DFl 


3.3 0.85 


3.3 0.83 




vdW-DFl(C09) 


3.0 1.06 


3.0 1.04 




vdW-DFl(optPBE) 


3.15 1.01 


3.15 0.99 




vdW-DFl(PW86r) 


3.15 1.10 


3.15 1.07 




vdW-DF2 


3.25 0.72 


3.25 0.69 




vdW-DF2(C09) 


3.1 0.65 


3.1 0.63 




Exp. 


0.85 6 






"Ref. 88 6 Rcf. 89 



tent with our results. Binding separations are unavail- 
able for any of the considered systems, complicating a 
clear assessment of the functionals. vdW-DFl typically 
overestimates binding separations by about 0.2-0.3 A for 
purely van der Waals bonded systems . 5 ' 7 ' 21 We find that 
the other variants reduce the separations, in particular 
those using C09 exchange. 

Caciuc and coworkers^ also calculated benzene ad- 
sorption on graphene and BN both for vdW-DFl and 
vdW-DF2. Our results agree well with theirs, with the 
exception their separations are about 0.1 A shorter for 
benzene on BN than ours. Using with semi-empirical 
DFT-D2j2i they found the same optimal sites as us. The 
fact that their study is based on a projector-augemented 
waves implementation in VASP^ while ours rely on ul- 
tasoft pseudopotential and DACAPO,^ 5 strengthens our 
confidence in the overall precision of both our and their 
study. The adsorption results for benzene on graphene 
also closely match the results in the seminal study by 
Chakarova-Kack and coworkers.— 



B. Induced dipole 

The upper panel of Fig. [S] visualizes the charge transfer 
density induced by the adsorption of C60 onto graphene, 
as given by Ap(r) = -eAn(r) = pc60-gr ap hene(r) - 
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h [A] 

FIG. 5: Upper panel: Charge transfer isosurfaces for C60 
on graphene Ap(r). Blue is negative, red is positive. Lower 
panel: Dipole induced as a function molecule-to-sheet separa- 
tion. Upper full black curve is for C60 on graphene, lower full 
blue is for C60 on boron nitride. Upper dashed is for benzene 
on graphene and lower is for benzene on BN. 



Pgraphcnc 

(r) — pc6o(r)' In the DFT calculations, grid sen- 
sitivity was minimized by using the same the same atomic 
coordinates relative to the underlying grid for the full sys- 
tems as for the separate fragments. The induced charge 
transfer density has a rich dipolar structure. The mag- 
nitude of this dipole depends on the separation between 
the molecule and the sheet. 

The lower panel of Fig. [5] shows the induced dipole 
as function of the benzene-to-sheet (dashed) and C60- 
to-sheet separation (full) for both graphene (black) and 
BN (blue) sheets. At typical binding separations and 
beyond, only C60-on-graphene produces a non-negligible 
dipole. The dipole that forms between C60 on BN is 
tiny, yet if the molecule is pushed towards the surface, 
the magnitude of the dipole increases; while for benzene 
adsorbants, the magnitude of the dipole barely changes 
with separation. 

To get a crude estimate for the total charge trans- 
fer from C60 to graphene, we first project the charge 
transfer density Ap(r) onto the z-axis Ap(z). We next 
identify the position zq where Ap(zo) = that lies be- 
tween the adsorbant and the sheet Ap(zo) = 0. Finally, 
we integrate over the projected charge transfer density 



smaller than z$: AQ — f*° dzAp(z). The estimated 
charge transfer, of about AQ « 0.04e, is quite insensitive 
to binding separation, just like the induced polarization 
^ind = J dz zAp(z) (as shown in the lower panel). 

The estimated charge transfer is moderate compared 
to 0.3 e found for tetrafluoro-tetracyanoquinodimethane 
Ci2(NF)4 on graphene^ This molecule, basically a 
hexagonal carbon ring with two arms of (NC)2 and four 
fluorine (F) in place of H, has been proposed as a p- 
dopant on graphene. That this systems results in a sig- 
nificant charge transfer is not surprising considering the 
electronegativity of fluorine. Since we find some charge 
transfer for C60 and none for benzene, we speculate that 
functionalized C60 might be a better candidate than 
functionalizcd benzene for doping graphene. 



IV. ANALYSIS OF EXCHANGE AND 
CORRELATION PART OF VDW-DF BINDING 

We analyze the role of exchange and correlation in the 
vdW-DF account of adsorption of benzene and C60 on 
graphene. This analysis benefits from comparing these 
two molecules, since they have a similar interface to the 
sheets, but different sizes. This property is first used to 
illustrate the inherit failure of local approximations to 
the exchange-correlation functional, as they are unable 
to capture the increased attraction of C60. Next, we 
examine the role of different exchange and correlation 
components to the potential energy curves of benzene 
and analyze the adsorption energy trends in light of these. 
We also confirm that the bulk of the difference between 
the PECs curves of benzene and C60 adsorption arise 
from non-local correlation. A more detailed analysis of 
the non-local correlation of vdW-DF follows including 
the role of density separations and density regimes. 



A. Failure of local approximations 

The fact that LDA does not include van de Waals 
forces is well established, yet one can still find recent 
LDA-based studies for sparse matter. Perhaps this is 
due to lack of awareness of DFT-D or non-local correla- 
tion functional methods. Perhaps a pragmatic point of 
view is taken: The underlying cause of the binding can 
be irrelevant for a given investigation, often the case for 
band structure calculations^ 

For some van der Waals bonded systems, an LDA ac- 
count of exchange-correlation may predict binding ener- 
gies and separations that compare well with experimen- 
tal. This occurs because the LDA exchange (as well as 
correlation) part is inadequate for sparse matter^ once 
replaced with a more sophisticated GGA functional, most 
or all this exchange binding vanish. We illustrate how lo- 
cal approximations to exchange-correlation functional is 
incapable of capturing the larger binding energy of C60 
compared to benzene on graphene. 
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FIG. 6: Potential energy curves for benzene (bz) and C60 on 
graphene for LDA, a modified LDA account, and vdW-DF. 




FIG. 7: Density isosurfaces for benzene on graphene at 3.0 A. 
The outer isosurface is for 0.006/afJ ohr and the inner (yellow) 
is for 0.027. They encapsulate all but 1.6% and 8% of the 
electrons. Figurge generated with VMD.— 



Figure [5] compares the potential energy curves for C60 
(black full curves) and benzene (orange full curves) . The 
full thick curves are obtained using vdW-DFl while the 
two upper thin curves are with LDA. It produces binding 
separations that are within the range of what the differ- 
ent version of vdW-DF predict. The binding energy is 
0.25 eV, about half of the experimental. Most revealing, 
however, is the puny 9% increase in LDA binding energy 
when comparing with C60 adsorption. 

As Fig. [5] confirms, LDA does not capture the long- 
range van der Waals forces. Neither would of course any 
other local approximation to the exchange-correlation 
functional. To illustrate how a fitted local approxima- 
tion might fail, we consider the following toy exchange- 



correlation functional, in the spirit of Xo£^, 

^ LDA = (l + a)£*c DA - (5) 

Here a is taken, for the sake of the argument, as a fitting 
parameter, which we set to 0.175 to produce the exper- 
imental adsorption energy of benzene on graphene, with 
corresponding PEC given by the dashed orange curve in 
Fig-El The optimal separation of about 3.0 A is likely an 
underestimation, but it is not that different from 3.25 A, 
as predicted by vdW-DFl(C09). However, when using 
this fitted description for a system of different size, the 
limitations, as also discussed in Ref. [H, become evident. 
Applied to C60 on graphene, the same toy functional pro- 
duces the black dashed curve of Fig. [5] The correspond- 
ing adsorption energy is merely 8% larger than that for 
benzene, completely off the experimental increase of 70%. 

Fig. [7] shows density isosurfaces for a benzene-to- 
graphene separation of 3.0 A. The well separated inner 
isosurfaces encapsulate 92 % of the total density density, 
while the outer, which is connected in some regions at the 
interface, encapsulate all but 1.6% of the total density. 
This small density overlap, even at such a short separa- 
tion, gives a clear indication why it is futile to capture 
long van der Waals interactions with local density func- 
tionals, in which the binding must be described alone by 
the local density magnitude. 

The presented LDA analysis also indicates the problem 
of combining accounts of van der Waals forces with GGA 
exchange functionals that induce some exchange binding, 
in particular if fitting to unrepresentative training sets 
are used. For bigger systems, in the interface-normal 
direction, the balance between repulsive and attractive 
terms can be shifted because a GGA exchange attraction 
only scales roughly with the size of the interface and it 
cannot keep up with the increased non-local attraction. 

We note that the popular PBE functiona l 49 ' 68 ' 69 has a 
certain unphysical exchange attraction and that it is of- 
ten used together with pair-potential accounts of van der 
Waals forces. The optPBE functional may also induce a 
small exchange attraction at moderate separations, and 
is in this sense not optimal. Formally, so does revPBE, 
but this exchange functional is so repulsive that this is 
not an issue in practice. 

B. Role of exchange and correlation version 

Figure [S] presents the exchange and correlation compo- 
nents of the PEC (upper panel) and their derivates (lower 
panel) for benzene on graphene as obtained with differ- 
ent versions of vdW-DF. The values presented earlier in 
table U are discussed in terms of these curves. The upper 
blue curve in the upper panel (circular dots) shows the 
revPBE exchange part of the vdW-DFl PEC, hereafter 
referred to as the revPBE curve; and in the same man- 
ner for the other exchange curves. The revPBE curve is 
above the other exchange curves. Hence, in agreement 
with table HI canonical vdW-DFl (with revPBE) gives 
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FIG. 8: Upper panel: Exchange and correlation energy com- 
ponents to the potential energy curve (PEC). Lower panel: 
Derivatives PEC components. 



the smallest binding energy. The revPBE curve is also 
the least steep among the exchange curves, as shown by 
its derivative in the lower panel. This result is related to 
revPBE producing the largest optimal separations. 

The C09 curve (green with square markers) is the 
steepest explaining why vdW-DFl(C09) and vdW- 
DF2(C09) produce the shortest separations. The C09 
curve becomes similar to the revPBE curve at 3.8 A. 
This similarity reflects the similar enhancement factors 
F x of revPBE and C09 in the large s-limit, as mentioned 
in section Hi A[ 

The optPBE curve has about the same depth as the 
C09 (upper panel), but is less steep (lower panel). Hence 
the adsorption energy of vdW-DFl(optPBE) is similar 
to that of vdW-DFl(C09), while the optimal separation 
is larger. 

The revPBE, optPBE, and C09 curves have a posi- 
tive value at a separation of 3.8 A. Since these curves 
inevitably will decline towards zero as the overlap be- 
tween the fragments vanish, these exchange variants have 
the potential to induce spurious binding at large separa- 



tions; a result related to the saturation of their enhance- 
ment factor F x (s) at large separations. 49 Whether spuri- 
ous binding indeed arise for a given binding curve depend 
on the balance between the different terms contributing 
to a given PEC. 

The PW86r curve stays below at all separations, thus 
this exchange choice does not induce binding for this sys- 
tem. For small molecular dimers, this non-binding prop- 
erty has been shown to be a consequence of the fact that 
the PW86r enhancement factor F x (s) always increases 
with s42 Our result indicates that this non-binding prop- 
erty carries over to to extended systems, which further 
motivates choosing this exchange partner. 

The full curves in the bottom part of the upper panel 
of Fig. [8] show the non-local correlation contributions to 
the PEC. The dashed include contributions from LDA 
correlation. 

The vdW-DFl correlation curve is significantly deeper 
than the one for vdW-DF2. This result relates to vdW- 
DFl's larger adsorption energies. The steepness of these 
two curves, shown in the lower panel, however differ much 
less from each other than that the steepness of the ex- 
change curves. Thus, within the selected range of func- 
tional variants considered here, the binding separation is 
primary set by which exchange variant is chosen. 

Table Q] shows that for a given exchange variant, vdW- 
DF1 predicts merely about 0.05 A shorter binding sepa- 
rations than vdW-DF2 despite the much larger binding 
energies. On the one hand, vdW-DFl and vdW-DF2 cor- 
relation greatly affects the binding energy, while on the 
other hand, the correlation version matters less for the 
binding separation. This observation makes the bind- 
ing separation a better parameter to assess the quality 
of an exchange partner to vdW-DF than the energy. In 
the same vein, if a strategy of fitting exchange to high- 
accuracy data sets is chosen, we expect the binding sep- 
aration to be the better target. 

The extremely steep LDA exchange curve (dotted 
curve in upper pane) seems to explain why LDA for 
some systems produce binding energies comparable to 
a van der Waals bond but for others severely underes- 
timate it. At larger separations, the curve goes rapidly 
to zero as the overlap between the fragments vanish. At 
short, the curve is even deeper than the vdW-DFl cor- 
relation curve. A slight shift in the kinetic-energy repul- 
sion or electrostatic interaction could move the molecule 
slightly inwards or outwards, adjusting the adsorption 
energy greatly. 



C. Benzene versus C60 adsorption 

Figure [5] confirms our assumption that the reason why 
C60 binds more and at shorter separations than benzene 
arise from the increased van der Waals interaction. The 
upper curves gives AEq 60 — AEq z on graphene. AEq 
has a mostly-repulsive contribution the PEC of C60 and 
benzene. The lower gives the same for the non-local cor- 
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relation part AEf> C60 - AEf< hz . Since the two lower 
curves are far deeper than the two upper, the increased 
binding of C60 over benzene on graphene arise primarily 
from non-local correlation. Similar results are found for 
benzene and C60 on BN. 

The figure further shows that the non-local correlation 
of vdW-DFl deepens the PEC more than that of vdW- 
DF2 does, while the small shift from Eq is similar for the 
two functional. This result is related to the bigger role 
of non-local correlations in vdW-DFl. That the bind- 
ing energies of the two functionals do not differ more is 
related to that revPBE exchange functional is more re- 
pulsive than PW86r. As we go from benzene to C60, the 
repulsive component do not change much. Therefore, the 
larger van der Waals component of vdW-DFl compared 
to vdW-DF2, reduces the binding separation more, and 
increases the binding energy more. 

The result that van der Waals forces is the primary 
cause of the increased binding of C60 compared to ben- 
zene on graphene (and BN) shows that these systems 
serve as a good starting point to explore properties these 
forces. Results for these systems can highlight the role 
played by different length scales or density separations in 
the non-local correlation. 



D. Role of different density separations 

To elucidate the role size and geometry have on the 
van der Waals attraction, we also study how much var- 
ious density separations contribute to the vdW-DF in- 
teraction energy. We perform this analysis for benzene 
and C60 on graphene and focus on the non-local correla- 
tion. For this purpose, we introduce a lower and upper 
separation cutoff in the vdW-DF kernel: 

fo,A(ri, ra) = 0(r x , r 2 )9[(R + A) - r 12 ]6[r 12 -(R- A)] , 

(6) 

where 8 is the Heaviside function. Just like for a full 
calculation, we subtract off, for a given kernel, the non- 
local correlation energy of the separate fragments (A,B), 

= *%R, a [Full] - A [A] - E* KA [B] . (7) 

This subtraction allows us to obtain a Separation- 
decomposed measure of the Non-Local correlation con- 
tribution to the Interaction energy (SNLI). This ap- 
proach is similar in nature to the non-local correlation en- 
ergy density analysis used by Lazic and coworkersi 75 ' 91 i 97 
However, while their analysis has emphasis on which spa- 
tial regions dominate the response, ours focus on the role 
of density separations . 

The columns of Fig, QJJ] show the SNLI for separation 
between R - A and R + A, with A = 0.1 A. The upper 
panel is for benzene on graphene and the lower is for C60 
(and benzene once more for comparison). The curves are 
guides to the eye. In all calculations, the adsorbant has 
an atomic separation of h — 3.6 A from the graphene 
sheet, the optimal one for benzene with vdW-DFl. At 



separations smaller than about 1 A, the SNLI is positive, 
while it is negative beyond. Its oscillatory shape reflects 
the shape of vdW-DF kernel, displayed in Fig. [5] The 
SNLI curves in Fig. [TU] are dominated by negative values 
because we have subtracted off the interaction within iso- 
lated fragments; short separations only contribute at the 
interface between the two fragments. 

The SNLI of vdW-DFl and vdW-DF2 are very similar 
for small density separations R in these systems. Be- 
yond the repulsive peak at 0.7 A, they start to deviate, 
but first at about 1.7 Ado this deviation become substan- 
tial. This result is linked to the stronger asymptote of 
vdW-DFl compared to vdW-DF24iiik 9 £ The similarity 
for small R reflects the fact that their description of the 
plasmon dispersion u; q [underlying the non-local correla- 
tion of expression (JJ] becomes similar for large gi 29 i 32 

Comparing the benzene (full) and C60 (dashed) curves 
in the lower panel, we see that they are almost identical 
for short separations: The benzene-graphene and C60- 
graphene interface is essentially the same as far as the 
non-local correlation is concerned. Beyond about 2 A, 
the benzene and C60 curves start to differ noticeably, 
and at 8 A, the contribution to the non-local correlation 
interaction is more than four times larger for C60 than 
for benzene, both in the case of vdW-DFl and vdW-DF2. 
Thus, if we gradually consider larger systems (compared 
to the size of interface), the difference between the non- 
local correlation of vdW-DFl and vdW-DF2 increases. 
This observation is clearly also related to the stronger 
asymptote of vdW-DFl. 

The separation analysis of Fig. QJj] also highlights the 
strikingly different nature of a vdW-DF and a DFT-D 
account of van der Waals forces. A similar histogram for 
DFT-D, or strictly pair-potential account for that mat- 
ter, would begin at the atomic separation, which in this 
case is 3.6 A. In contrast, and even despite the extended 
nature of the systems considered, most of the non-local 
correlation interaction energy in vdW-DF is already ac- 
counted for at separation shorter than 3.6 A. It is not 
surprising that vdW-DF is enhanced at moderate sep- 
arations compared to a corresponding asymptotic pair 
potential account as discussed in Refs.[H|0- This enhance- 
ment is a reflection of the plasmon- nature of the vdW-DF 
design. 



E. Analysis of density regimes 

Our analysis shows that most of the van der Waals in- 
teraction in vdW-DF arise from separations shorter than 
the ionic separation between the molecules. Because low- 
density regions are closer to each other, this result sug- 
gests that the low-density regions contribute the most 
to the interaction energy arising from the non-local cor- 
relation energy of vdW-DF. To quantify the role of low- 
density regions, we introduce a low-density cutoff in eval- 
uating the non-local correlation of Eq. (JJJ) (rather than 
separation cutoffs). By gradually increasing this cutoff, 
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FIG. 9: The difference between the components of the PEC 
for benzene and C60 on graphene, shown for vdW-DFl (blue, 
full)) and vdW-DF2 (red,dashed). 
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FIG. 10: vdW-DF non-local correlation contributions to in- 
teraction energy as a function of separation between densities 
(SNLI) as specified by Eq. (0. The transparent histograms 
show the result for benzene (bz) on graphene in the upper 
panel for respectively vdW-DFl (blue) and vdW-DF2 and 
C60 on graphene in the lower. In both curves the full curves 
indicate the result for bz on graphene, while in the lower the 
dashed curves indicate the result for C60 on graphene. 
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FIG. 11: Fraction of non-local correlation interaction that 
remains as a low-density cutoff is introduced, expressed as a 
function of the fraction electrons removed. The blue (red) 
square dots is for benzene at a molecule-to-sheet separation 
of ft = 3.0 A using vdW-DFl (vdW-DF2), while the circular 
dots indicates the same for a separation of h = 3.6 A. The 
insert relates the low-density cutoff A cut (in terms of atomic 
units) to the fraction of density removed. 



we systematically gauge the role of various-density re- 
gions. 

Figure [IT] shows the result of this density analysis for 
benzene on graphene at two representative separation of 
h = 3.0 A and h = 3.6 A. The vertical axis indicates 
the fraction of the non-local correlation interaction en- 
ergy AZ?" 1 that remains after a given low-density cutoff 
has been introduced. The same density cutoff is used in 
the full calculation as in the reference calculations. The 
horizontal axis indicates the fraction of the total num- 
ber of electrons that has been removed in the full system 
for a given density cutoff. The insert relates the density 
cutoff [in atomic units] (vertical axis) to the fraction of 
electrons removed (horizontal axis). 

The figure reveals that the vdW-DF interaction energy 
is very sensitive to low-density regions. For the case of 
h =3.0 A , removing 2% of the density (as indicated in the 
outer contour in Fig. [7]) removes 15% of the interaction 
energy; removing 8 % of the density (indicated in the 
inner counter of Fig. [7]) removes 65% of the interaction 
energy. If we remove 50 % of the density, starting from 
the smallest density, 97% of the interaction energy is re- 
moved. For a vdW-DF like functional that weighted den- 
sity equally, the non-local correlation would be quadratic 
in the density and the corresponding numbers would be 
4%, 15%, and 75%. These widely different trends reflects 
how the kernel (j>[n) itself is a functional of the density n. 

Low density regions are even more important at the 
larger separation h = 3.6 A. To crudely explain this 
effect, we can consider that when the molecules move 
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further apart, the density between the molecules is low- 
ered, because the density overlap of the molecules is re- 
duced. At the larger separations, we also find that vdW- 
DF1 is somewhat more sensitive to low density regions 
than vdW-DF2, while at 3.0 A the curves are very sim- 
ilar. This trend agrees well with the separation analysis 
shown Fig.[TUl since vdW-DFl and vdW-DF2 differ more 
at larger separations. 

Our analysis show that interactions between different 
low density regions (as connected by the two-point kernel 
<p(r,r')) dominate the contributions to the interaction 
energy. Clearly the response of low-density regions is 
much stronger that that of high density regions. Since 
part of the low-density regions are closer to each other 
(the region close to the interface) this result can be linked 
to the importance of shorter density separations (than 
atomic separations between fragments) established in the 
previous subsection. 

Finally, we note that the importance of low-density 
regions has an important bearing on how we build our 
intuition of how strong a van der Waals interaction will 
be. According to the presented analysis, it is not the 
number of atoms that roughly scale with strength of 
van der Waals interactions but more the amount of low- 
density regions in a fragment. This observation can be 
for instance linked to the offset in the binding trends of 
PAHs— and alkanes^i on graphene and the big re- 
sponse of hydrogen atomsi 5 ' 100 



V. CORRUGATION RESULTS 

The predicted corrugation is very sensitive to which 
version of vdW-DF and exchange variant is used. Cor- 
rugation is here understood as the potential energy vari- 
ation of a molecule on different sites on the surface. We 
have recently shownii for H2 on Cu(lll) that the cor- 
rugation is even more sensitive to switching between a 
vdW-DF to a DFT-D account. The adsorption energies 
and optimal separations, presented in section Hm is also 
significantly affected by the vdW-DF variant used; for 
benzene on graphene, the optimal separations differ by 
8% and the energy by 70%. Nevertheless, the vdW-DF 
sensitivity to corrugation stands out: for this system, the 
diffusion barriers can change by up to a factor of three as 
we switch the exchange variant and correlation version. 

This huge sensitivity arise because the corrugation de- 
pends exponentially on the binding separation; slight 
changes in separation greatly affect the corrugation. 

Figure Q2] shows the how the optimal binding energies 
varies for the different high-symmetry sites on graphene 
and BN for benzene and C60 adsorbants, as calculated 
with the considerent variants of vdW-DF. It shows that 
the corrugation depends strongly on the correlation ver- 
sion used and in particular on the exchange choice. The 
upper panels show the results for benzene on graphene 
and on BN. The lower, show those for C60 with the car- 
bon hexagon facing graphene and BN, which is preferred 



over the pentagon face. Graphene sites are identified in 
Fig. [3J the BN sites are similar, but with twice as many 
top sites: the hexagon can be centered on N or on B 
(top-N and top-B). The vertical scale in the lower pan- 
els are twice those of the upper; the corrugation of C60 
adsorbants is roughly twice that of benzene. 

Benzene can easily diffuse on graphene since the 
molecule can pass between top and bridge sites which 
have very similar energies. This open-path dynamics is 
similar to that of benzene on Cu(lll) i 18 i 47 On BN this 
path is closed, since from a top-N site, the molecule must 
either move to a top-B site or a hollow (or between) to 
get to the next top-N site. The corrugation energies are 
also generally somewhat larger. Thus our calculations 
indicate that benzene will freeze at much higher tem- 
peratures on BN than on graphene. C60 has a similar, 
though less pronounced, low-energy path on graphene. 
Rotational barriers are in general small both for benzene 
and C60. We note that since C60 adsorbants induce 
a dipole in particular on graphene, in-plane dynamics 
must be accompanied by continuous charge redistribu- 
tion. Since electron-hole excitations express such charge- 
relocations, this will open a damping mechanism causing 
frictioni£i"i2ii for the dynamics for C60 on graphene. 

There is a question of whether C60 might diffuse like a 
soccer ball on graphene. According to classical molecular 
dynamics calculations, this is not the case for thermally 
induced motionjA24 However, such a classical analysis fail 
to account for dynamics of the friction that we argue will 
be present. 

Comparing the functionals, we find that vdW-DFl pre- 
dicts the smallest corrugations followed by vdW-DF2, 
while vdW-DFl(C09) clearly predicts the largest. In 
between are vdW-DF2(C09), vdW-DFl(PW86r), and 
vdW-DFl(optPBE) which are more or less tied. The 
succession follows closely that of the binding separations 
from smallest to largest as listed in table |U For ben- 
zene, vdW-DFl clearly predicts a larger corrugation than 
vdW-DF2, but for C60 the difference is smaller. This 
trend is reflected in the more similar separations pre- 
dicted for C60 adsorption. 

The observation of a strong link between optimal sepa- 
ration and corrugation is somewhat clouded by fact that 
the separation is optimized separately for each site. To 
make a clearer analysis, we consider how the difference 
between the PECs of two specific sites (APEC) changes 
if we change the exchange variant or correlation version. 

Fig. Q2] shows, for benzene on graphene, that APEC 
barely changes with the exchange variant and correlation 
version used. In contrast, the full PEC might change 
much. Thus, the corrugation of graphene and BN is pre- 
dominantly set by the binding separation (given by the 
minimum of the full PEC), making the huge sensitivity 
to functional version arise indirectly as consequence of 
the exponential sensitivity to binding separation. 

In Fig. [T31 the indirect sensitivity to corrugation is doc- 
umented separately for the mostly-repulsive AEq part of 
the PEC (upper panel) and for the non-local correlation 
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FIG. 12: The variation in adsorption energy as a function of the site, for benzene (upper row) and C60 (lower row) adsorbants 
on graphene (left column) and on boron nitride (right column). The data points (connected by lines) are obtained with different 
variants of vdW-DF. 



part A££' (lower panel). In the upper panel, the differ- 
ence between AEq of the top site and that of the bridge 
and hollow, for the exchange choice of revPBE, PW86r, 
and C09. And in the same manner in the lower panel 
for AEf for vdW-DFl and vdW-DF2 correlation. The 
curves conform closely to each other (the same goes for 
optPBE exchange which is not shown). The small dif- 
ference in the AEo curve is much bigger that the dif- 
ference between the non-local correlation contribution 
curves AEf. The vertical axis of the lower panel spans 
one-tenth of that of the upper panel. 

The similarity of the tiny contributions to corrugation 
arising directly from vdW-DFl and vdW-DF2 correlation 
can be explained by the similarity of their SNLI curves 
shown in Fig. [10] for short separations. The curves in 
the upper panel depend essentially exponentially on sep- 
aration — for the largest separations noise contributes 
to the tiny difference between the top and hollow PEC. 
The curves document how corrugation is set by the opti- 
mal separation. It is the effect the exchange-correlation 
variant has on the binding separation that causes the 



corrugation to be strongly affected by the variant. 



VI. TOWARDS A SURFACE-BASED 
EXPERIMENTAL REFERENCE DATABASE 

Adsorption on clean surfaces provide outstanding op- 
portunities to compare theory and experiments. Spe- 
cific molecules can be carefully deposited on surfaces 1 ^ 
and imaged with scanning-tunneling (STM) or atomic- 
force microscopy (AFM). Several experimental probes 
can be used to gain detailed quantitative insight that 
can be compared with calculated numbers. One exam- 
ple is the comparison of theoretical and experimental 
vacuum-level shifts that can be used to deduce binding 
separations and charge transfer of weakly chemisorbed 
molecules J£§rl22. Another is the measurements of reso- 
nance levels arising in back-scattering of light molecules, 
which provide a powerful connection between experi- 
ments and theory . 46 ' 47 i 109 ~ 21£ These experiments can be 
used to deduce the shape of the PECs as well as corruga- 
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FIG. 13: Corrugation sensitivity to separation h between ben- 
zene adsorbant and graphene sheet. Upper panel: Difference 
between AEo(h) for a given site and the top site for three dif- 
ferent exchange functionals. Lower panel: In the same man- 
ner for the non-local correlation part A_E°' f° r vdW-DFl and 
vdW-DF2. 



tion and could give detailed information on interactions 
at many different length scales. 



A. Proposed concerted theory-experiment 
approach 

Based on our results and analysis, we propose a con- 
certed theory-experiment approach to probe the magni- 
tude of the various exchange-correlation contributions 
to physisorption on surfaces. The proposed approach 
is complimentary to others that emphasize a back- 
scattering based compariso n 46 ! 47 ' 109- -^ (limited to light 
molecules) or vacuum-level shift comparison^^r-i ^ (lim- 
ited to systems with a sizeable charge transfer at the 
interface). 

The suggested approach takes three steps: First, one 
compares experimental and calculated corrugation and 
thus determine the binding separation. Second, based on 
the separation, one estimates how much of the adsorp- 
tion energy arise directly from the exchange-correlation 
energy. Third, comparing the adsorption of molecules 
with the same interface geometry but different shapes 



and volume, like benzene and C60, one probes how much 
of the adsorption energy arise from the van der Waals 
interaction. 

The first step is motivated by the fact that the pre- 
dicted corrugations are strongly sensitive to the vdW- 
DF variant — in particular to the exchange (Fig. IT21) 
while at the same time the difference between PECs of 
two sites is very insensitive to the variant used fFig. [T5|) . 
This result indicates that the strong corrugation sensi- 
tivity, at least for benzene and C60 adsorption on BN 
and graphene, arise almost exclusively from the shift 
in separation induced by the replacement of exchange- 
correlation functional. Turning this around, we can can 
conclude that the corrugation essentially determines the 
separation. If the corrugation can be measured, DFT cal- 
culations can determine the binding separation, since the 
specific exchange-correlation used for this comparison is 
unimportant. 

Many other systems might also exhibit this kind of di- 
rect exchange-correlation insensitivity. The identification 
of additional systems of this kind merits further study, 
in particular if the corrugation has been accurately mea- 
sured. The validity of the argument presented here is also 
limited by the fact that the exchange-correlation account 
affects the Kohn-Sham (KS) orbitals and thereby kinetic 
energy contributions. However, for purely van der Waals 
bonded systems, the direct part of this effect is likely 
small. 

The second step also builds on the assumption that the 
Kohn-Sham orbitals are not much affected by the specific 
exchange-correlation functional used. This assumption 
is appropriate for weakly-bonded systems^ If so, the 
total exchange-correlation contribution to the adsorption 
energy can be assessed if both the binding separation and 
energy is known. The separation also tells us how the 
exchange-correlation functional affects the steepness of 
the PEC, since the van der Waals forces must counter the 
kinetic-energy repulsion to make the molecule physisorb 
at a given separation. 

The third step builds on the result that for molecules 
with similar interfaces to the surface (or sheet) like ben- 
zene and C60, but with different geometries and size, we 
can single out the effect of van der Waals forces (Fig. \§§ . 
Such results give us insight on the delicate balance be- 
tween the exchange and correlation contributions to the 
binding. These results can be combined with results 
for the asymptotic interactions, since advanced meth- 
ods can compute these to high accuracy 37 ' 113 " -^ Such 
analysis tells us how the van der Waals forces are modu- 
lated at shorter separations compared to the asymptotic 

formf 5i40^5? 

How to accurately determine the corrugation experi- 
mentally is a question beyond the scope of this paper, 
but we note that the advancements in STM and AFM 
is astonishing. One might drag molecules across the sur- 
face and measure the potential energy landscape ; 119 ' 120 
or observe the diffusion of molecules on the surface i 121 ' 122 
Based on diffusion barriers, corrugation can be esti- 
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mated. 



B. Surface-based database 

Many adsorption systems may have the structural re- 
lation as benzene and C60 on graphene and BN and 
that could make them good discriminators of exchange- 
correlation accounts. Other experimental-systems, like 
scattering of light molecules on noble metal surfaces or 
systems with vacuum level shifts also provide excellent 
links between theory and experiment. As experimental 
methods continue to develop and new experimental data 
becomes available, these could be linked to DFT calcula- 
tions or parameterized in such a manner that it becomes 
straightforward to assess the performance of sparse mat- 
ter methods. The construction of PEC based on reso- 
nance levels of H 2 backscatterin g 109 ' 111 makes it simple 
to assess the theoretical methodsi 46 ' 47 The relation be- 
tween corrugation and binding separation discussed in 
section [V] provides a different example. 

If we identify and organize several of these theory- 
experiment links for surfaces we can build up a database 
that could provide invaluable for the continued advance- 
ment of sparse-matter methods. Today datasets based on 
high-accuracy quantum-chemistry calculations play such 
a role? 41 - - — but surface electrons have a different collec- 
tive nature than that of small molecular dimersi 30 ' 123 
Surfaces are also especially suited to extract a range of 
experimental properties. The corrugation sensitivity to 
binding separation is central to the proposed approach. 
Other quantities sensitive to geometry, like the HOMO- 
LUMO gap of adsorbed molecules , 124 ' 125 may help un- 
cover other links between theory and experiment useful 
for assessing sparse-matter methods. 

VII. SUMMARY 

In this paper, we identify a comparison of the adsorp- 
tion of benzene and C60 on inert surfaces like graphene 
and BN as a good arena to evaluate exchange-correlation 
accounts of van der Waals forces and to gain insight 
into the nature of physisorption. Our study of ben- 
zene and C60 on graphene have shown that vdW-DFl 
and vdW-DF2 accurately capture the increased van der 
Waals binding of the larger molecule. 



A detailed analysis of adsorption results, exchange- 
correlation components to the PECs, and corrugation 
have been performed. Based on this analysis, we have 
made a number of observations. We have shown that 
C60 induces a moderate dipole on graphene, while for 
the other systems dipoles are tiny. We have shown that 
the binding separation is more sensitive to the exchange 
variant than the correlation version of vdW-DF. This re- 
sult making the separation a better target for assessing 
or optimizing an exchange account. We have shown that 
the difference between benzene and C60 binding arise 
predominantly from the increased non-local correlation, 
displaying fundamental shortcomings of relying on un- 
physical exchange binding. With our density separation 
decomposition of the non-local correlation interaction, we 
documented that vdW-DFl and vdW-DF2 are similar for 
short density separations, but differ much at large. The 
similar density-cutoff analysis reveals how the non-local 
correlation of vdW-DF is very sensitive to low-density 
regions. 

Further, we have demonstrated that the corrugation 
is very sensitive to the vdW-DF variant used, but also 
that a particular low-energy path exists for benzene on 
graphene, making it easy for this molecule to diffuse. The 
corrugation sensitivity have been shown to arise almost 
exclusively as an indirect effect arising from the shift in 
binding separation caused by the exchange choice and 
correlation versions. This result has lead us to suggest 
a experiment-theory approach to assess variants of, and 
gain insight into, the exchange-correlation functional and 
its account of van der Waals forces. 

The novel analysis approaches presented here, like 
density-separation and density-regime decomposition, 
can also be applied to other material systems and func- 
tional. Such approaches can allow us to view materials 
and methods in new ways helping the development of 
non-local correlation functionals 
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